function [fftdata, f] = fast_fourier_transform(data, Fs)
% 幅值频谱计算
%---input---
%data：数据，size1采样点，size2通道
%Fs：采样率
%---output---
%fftdata: 单边频谱，行为采样点，列为通道
%f: 频率横轴

T = 1/Fs;
L = size(data,1);
Y = fft(data);
P2 = abs(Y/L);
P1 = P2(1:L/2+1,:);
P1(2:end-1,:) = 2 * P1(2:end-1,:);
f = Fs*(0:(L/2))/L;
fftdata = P1;